Magnetoelectric nature of skyrmions in a chiral magnetic insulator Cu20Se03 
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Dielectric properties were investigated under various magnitudes and directions of magnetic field 
(H) for a chiral magnetic insulator Cu20Se03- We found that the skyrmion crystal induces electric 
polarization (P) along either in-plane or out-of-plane direction of the spin vortices depending on the 
applied if-direction. The observed if-dependence of P in ferrimagnetic, helimagnetic, and skyrmion 
crystal state can be consistently described by the d-p hybridization model, highlighting an important 
role of relativistic spin-orbit interaction in the magnetoelectric coupling in Cu20Se03- Our analysis 
suggests that each skyrmion particle can locally carry electric dipole or quadrupole, which implies 
that the dynamics of skyrmions are controllable by the external electric field. 

PACS numbers: 75.85.+t, 75.70.Kw, 77.22.-d 



Skyrmion is a vortex-like spin swirling object with 
particle nature [lml, which has recently been identified 
in chiral magnets 0, In metallic system, skyrmions 
can respond to the electric current through spin transfer 
torques from conduction electrons [6|-l8| ; Such electric con- 
trollability, as well as its partcile-like nature and nano- 
metric size, highlights skyrmion as a promising building 
block for next generation of magnetic storage device 0|. 
In chiral magnets, magnetic interaction acting on two 
neighboring spins Si and Sj consists of two terms; Si • Sj- 
like exchange interaction and Si x 5j-like Dzyaloshinskii- 
Moriya (D-M) interaction. Skyrmions are stabilized in 
the limited window of temperature (T) and magnetic 
field (H) over the helimagnetic ground state, and can 
appear as independent particle or in hexagonal lattice 
form within a plane normal to H (Fig. 1(c)) [J, While 
the previous observation of skyrmions has been limited 
to specific metallic alloys with chiral B20 structure such 
as MnSi0, FeGepjj, and Fei_ a; Co a .Si@,Ii3, recently the 
formation of skyrmion crystal (SkX) has also been iden- 
tified in an insulator Cu20Se03 through Lorentz trans- 
mission electron microscopy experiments on thin film[l2| 
as well as small angle neutron scattering study on bulk 
crystal [H Q . 

The crystal structure of Cu2 0Se03 belongs to the chi- 
ral cubic space group P2 X 3 (Fig. 1 (a)) fl5l4l7| . The 
magnetic ground state is proper screw, where spins ro- 
tate within a plane normal to the magnetic modula- 
tion vector q (Fig. 1(d)). In the bulk Cu 2 OSe0 3 , the 
SkX state is stabilized in the narrow T- and iJ-region 
called 'A-_phase' just below magnetic ordering temper- 
ature T c [l2|-[l3- The magnetic modulation period in 
the helical and SkX state (~ 630 A[H, 13 ) is much 
longer than the crystallographic lattice constant (~ 8.9 
A [la fl7j). and three- up one-down type local ferrimag- 
netic arrangement between two inequivalent Cu 2+ sites 
has been proposed [l7l4l9| . Notably, Cu20Se03 hosts 



magnetically-induced electric polarization (P) in the fer- 
rimagnetic, helimagnetic, and SkX state[12|. The cou- 
pling between helical spin texture and ferroelectricity 
has recently been reported for several compounds like 
TbMnO 3 [20|-[23|, which often enables unique magneto- 
electric (ME) response such as a control of magnetism by 
external electric field (E) . Although such a ME 

control has never been applied to particle-like spin object 
apart from the current-drive of spin vortex motion in a 
microdisk magnet, the reported emergence of P in the 
SkX state implies the possible manipulation of skyrmions 
by external electric fields in insulators. At this stage, the 
microscopic origin of ME coupling in Cu20Se03 has not 
been identified |l7l l27| , and the dielectric properties in the 
SkX state were investigated only for H || [111] [12]. To 
fully address the possibility of i^-control of skyrmions in 
insulators, further study of their ME properties is highly 
desirable. 

In this Letter, we investigated the dielectric properties 
of bulk Cu20Se03 single crystal under various magni- 
tudes and directions of H. The observed development 
of P in all magnetic phases can be well reproduced by 
the recently proposed d-p hybridization mechanism [28- 
[30|, suggesting an important role of spin-orbit interac- 
tion on the ME coupling in Cu20Se03. In the crystal- 
lized form, skyrmions are found to induce P along the in- 
plane or out-of-plane direction of spin vortices depending 
on the i^-direction with respect to the crystallographic 
axes. Our analysis shows that each skyrmion particle lo- 
cally carries electric dipole or quadrupole, which suggests 
that the dynamics of an individual skyrmion particle is 
controllable by the external electric field in insulators. 

Single crystals of Cu20Se03 were grown by the chem- 
ical vapor transport method [27[ • They were cut into a 
cubic shape with faces parallel to (110), (111) or (001), on 
which silver paste was painted as electrodes. To deduce 
P, we measured the polarization current at constant rates 
of i7-sweep (10.6 Oe/sec) or ^-rotation (l°/sec), and 
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FIG. 1: (Color Online), (a) Crystal structure of Cu 2 OSe0 3 
with two distinct Cu 2+ sites with different oxygen coordi- 
nation, (b) Symmetry elements compatible with the crystal 
lattice of Cu20Se03. (c) Magnetic skyrmion crystal formed 
within a plane normal to applied magnetic field (H). Back- 
ground color represents the out-of-plane component of local 
magnetization vector (m z ). (d) Proper screw helical spin tex- 
ture with a magnetic modulation vector q. In (c) and (d), 
symmetry elements compatible with each spin texture are also 
indicated; Green dashed lines represent two- fold rotation axes 
(2) or two- fold rotation axes followed by time-reversal (2 r ), 
and a small black hexagon does a six-fold rotation axis (6) 
along the out-of-plane direction, (e) Schematic illustration of 
d-p hybridization mechanism, (f)-(h) Magnetically-induced 
electric polarization (P) under various directions of H for 
Cu20Se03, predicted by the symmetry analysis (see text). 



integrated it over time. Here, polarization current was 
measured with an electrometer without applied E. Since 
the transition from paraelectric to ferroelectric state in 
the process of zero-field-cooling (ZFC) should produce 
equal population of P-domains with no net polarization, 
we assumes P = for the initial H = state just af- 
ter the cooling without E and H. Magnetization M and 
ac magnetic susceptibility x' ( a ^ 700 Hz) were measured 
with a SQUID magnetometer. 

To investigate the microscopic origin of ME coupling 
in Cu2 0Se03, we first measured the indirection depen- 
dence of P in the collinear (i.e. ferrimagnetic) spin state 
at 2 K with H = 0.5 T. Figures 2(a) and (b) indicate 
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FIG. 2: (Color Online), (a) [110] and (b) [001] components 
of P simultaneously measured in H rotating around the [110] 
axis, (c) [001] component of P under H rotating around the 
[001] axis. Both measurements were performed at 2 K with 
H — 0.5 T (i.e. collinear ferrimagnetic state). Dashed lines 
indicate the theoretically expected behaviors based on Eq. 
(1), and arrows denote the direction of ii-rotation. In (d) 
and (e), the experimentally obtained relationships between 
the directions of P and M in the ferrimagnetic state, as well 
as the definition of 6 and <f> (the angle between the P-direction 
and the specific crystal axis), are summarized. Here, the di- 
rections of M and P are indicated, while x denotes P = 0. 



the [110] and [001] component of P (P[no] and P[ooi]) as 
a function of indirection, simultaneously measured with 
two sets of electrodes. Here, H rotates around the [110]- 
axis, and 6 is defined as an angle between indirection 
and the [001] axis (Fig. 2(d)). We also measured the de- 
velopment of P[ooi] f° r H rotating around the [001] axis 
(Fig. 2(c)), where an angle between H and the [100] axis 
is defined as <\> (Fig. 2(e)). As a function of respective 
iY-rotation angle, both P-profiles show sinusoidal modu- 
lation with the period of 180°. 

Recently, Jia et al. suggested that at least three mi- 
croscopic mechanisms can be considered as the sources 
of ME coupling 
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[29|. Two of them originate from 
the correlation between two neighboring magnetic mo- 
ments (mi) and (raj) ((•••) denoting the expected value 
in the magnetically ordered state), and they provides lo- 
cal electric polarization proportional to (rhi) • (rhj) and 
(rhi) x (mj)(31j, respectively. However, these two mech- 
anisms predict P = for the collinear spin state in 
Cu20Se03, due to the cubic symmetry of crystal lat- 
tice or to the relationship (rhi) x (rhj) = 0. The third 
mechanism of ME coupling, namely the d-p hybridization 
model [28l-[30[ l32| , arises from the interaction between a 
ligand (oxygen) ion and a transition metal (copper) ion 
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FIG. 3: (Color Online). Magnetic field dependence of (a) magnetization M, (b) ac magnetic susceptibility x \ an d (c) electric 
polarization P for Cu20Se03 measured at 57 K with H || [001]. The corresponding profiles for H || [110] and H || [111] are 
also indicated in (d)-(f) and (g)-(i), respectively. Letter symbols F, S, H, and H' stand for ferrimagnetic, SkX, helimagnetic 
(single ^-domain) and helimagnetic (multiple ^-domains) states, respectively. Dashed lines indicate the theoretically expected 
behavior based on Eq. (1) for the single domain helical spin state, (j) Development of helical spin structure under applied H. 
(k) P-M correspondence for the helical and SkX states in Cu2 0Se03, calculated based on the d-p hybridization model (see 
text). 



with a single magnetic moment (fhi) (Fig. 1(e)). This 
model predicts the local electric dipole pij in the form 
of oc (eij • (rhi}) 2 eij [30|, where is a unit vector 
along the bond direction[33j. In this scheme, the cova- 
lency between metal d and ligand p orbitals is modu- 
lated depending on the local magnetization direction via 
spin-orbit interaction, and thus the local electric dipole 
is produced along the bond direction [3J. For the given 
coarse-grained spin texture m(r), the d-p hybridization 
scheme approximately predicts P and M in the form of 

P oc -t^ J ' rh{r)) 2 eijdr (1) 

^ i,3 

and M oc j ra(r)dr/ j dr, respectively. Here, summation 
is taken over 80 Cu-0 bonds within a crystallographic 
unit cell and integral is over a magnetic unit cell. By 
assuming the collinear spin state with spatially uniform 
m(r), Eq. © predicts (P [110 ] , P [00 i]) = (2_Psin20,P(l - 
cos 20)) and P[ooi] = 2Psin2(/>, where P is a common 
amplitude of P. These formulas (dashed lines in Figs. 
2(a)-(c)) well reproduce the experimentally observed in- 
direction dependence of P, which strongly suggests the 
validity of d-p hybridization mechanism as the origin of 
ME coupling in Cu20Se03. The obtained relationship 
between the directions of P and M in the collinear spin 



state is summarized in Figs. 2(d) and 2(e). We can see 
P = for M || [001], while P || [001] for M || [110] and 
P || [111] for M || [111]. 

Next, we measured M and \' as functions of the mag- 
nitude of H at 57 K just below T c - 58 K, for H \\ [001] 
(Figs. 3(a) and 3(b)), H \\ [110] (Figs. 3(d) and 3(e)) 
and H \\ [111] (Figs. 3(g) and 3(h)). In the ground state 
(H = 0), proper screw spin order is realized with muti- 
ple g-domains with equivalent q || (001) directions due to 
high symmetry of cubic lattice [12Ml4| . Application of H 
first induces the alignment of q along H direction while 
keeping the screw-like spin structure, as detected as an 
enhancement of x /-va lu e around H = 100 Oe. Further 
increase of H induces continuous transformation of spin 
texture from proper screw to conical, and finally to col- 
llinear (Fig. 3(j)). It causes almost-linear increase and 
saturation of M- value. In the intermediate field region 
with 200 < H < 300 Oe, x' shows clear dip anomaly that 
signals the formation of the SkX state [12l4l4j. 

Figures 3 (c), (f) and (i) indicate the corresponding de- 
velopment of P under various directions of H. We found 
that all of ferrimagnetic, helimagnetic, and SkX states 
induces P || [111] for H || [111], P || [001] for H || [110], 
and P = for H || [001]. Such a relationship between 
P and H can be understood from the viewpoint of sym- 



4 



<g)H 




+1 



m z , p*, p 



■ 



>H|| [001] 
P = 

r-KllO] 

*[110] 

>H|| [110] 

i— [110] 
; [001] 

» H || [111] 

|— [110] 
*[112] 



|(c) : 



(e) 



(f) 



(g) 



FIG. 4: (Color Online). Calculated spatial distribution of 

(a) local magnetization vector m, (b)-(d) local electric po- 
larization vector p, and (e)-(g) local electric charge p for the 
skyrmion crystal state (see text). Magnetic field is along the 
out-of-plane direction, and the results for H || [001] ((b),(e)), 
H || [110] ((c), (f)), and H || [111] ((d), (g)) are indicated. The 
background color represents relative value of m z for (a), p z for 

(b) -(d), and p for (e)-(g), respectively. Here, m z and p z stand 
for the out-of-plane component of m and p, respectively. The 
dashed hexagon indicates a magnetic unit cell of the skyrmion 
crystal, or a single skyrmion particle. 



metry. Whereas the original crystal lattice of Cii20Se03 
belongs to non-polar space group P2i3 (Fig. 1(b)), the 
combination with the ferrimagnetic, proper screw (Fig. 
1(d)), or SkX (Fig. 1(c)) spin texture leads to symmetry 
reduction and may produce an electrically polar state de- 
pending on the indirection. In case of H || [111], only a 
three-fold rotation or screw axis along H remains unbro- 
ken and hence P \\ H || [111] may emerge (Fig. 1(h)) (l2| . 
For H || [110], only a screw (or screw with time-reversal) 
axis along the [001] direction survives and hence P || [001] 



is allowed (Fig. 1(g)). With H || [001], the orthogonal 
arrangement of screw (or screw with time-reversal) axes 
remains, therefore no P is expected (Fig. 1(f)). Such a 
situation is realized for all of ferrimagnetic (with H || M), 
proper screw (with H || g), and SkX states. The above 
predictions by the symmetry analysis are fully consistent 
with the experimental observation. 

Based on Eq. (pQ), we have further calculated the P- 
M profile for various spin textures in Cu20Se03. In the 
helimagnetic phase, H induces continuous deformation 
of spin texture from helical to collinear as shown in Fig. 
3 (j). The spin texture in this process is described by 

rn(r) oc e z M\ + [e x cos(<f • r ) + e. y sin(<f • r)], (2) 

where e x , e y and e z are unit vectors orthogonal to each 
other. e z and q are parallel to the applied indirection. 
Likewise, the SkX phase (Fig. 1(c)) is approximately 
given by [4] 

3 

m(r) oc 4M2+y^[4cos(g Q -f+7r)+4sin(g Q -f+7r)]. (3) 

a=l 

q a denotes one of three magnetic modulation vectors nor- 
mal to which forms an angle of 120° with respect to 
each other. e a is a unit vector orthogonal to e z and <f a , 
defined so that all q a • (e z x e a ) may have the same sign. 
Here, M\ and M2 scale with the relative magnitude of 
net magnetization along the indirection, and M\ — >• 00 
and M2 — >• 00 correspond to the collinear spin state; 
The prefactors for Eq. (j2j) and (j3j) are determined so 
as to keep the averaged spin density (j \ fh(r)\df/ f dr) 
fixed. When we define the P and M in the collinear 
spin state as Pq and Mq, Eq. (pQ) predicts the relation- 
ship P/Pq — |(M/M ) 2 - \ for the helimagnetic phase 
(Eq. (J2J)) as shown in Fig. 3 (k). In Figs. 3(f) and 
3(i), the experimentally observed P-H profiles are fitted 
with the combination of the calculated P-M correspon- 
dence (Fig. 3 (k)) and the measured M-H profiles (Figs. 
3(d) and 3(g)). The fitting curves (dashed lines) well re- 
produce the observed P-H profiles in the helimagnetic 
phase, including the sign change of P. P remains zero 
in the low-ii" region with multiple g-domains, probably 
due to the cancelation of P averaged over different do- 
mains. Likewise, we also calculated the P-M correspon- 
dence for the SkX phase with Eq. ^ (Fig. 3(k)). Our 
calculation predicts P = 0.41P for M = 0.46M , which 
roughly agrees with the experimentally observed values; 
P = 0.32P for H || [110] (Fig. 3(f)) and P = 0.33P 
for H || [111] (Fig. 3(i)) for the corresponding M- value 
in the SkX state [35|. The good agreement between the 
calculated and observed P-profiles, including the repro- 
duction of complicated sign change of P, confirms that 
the d-p hybridization mechanism adopted here as the mi- 
croscopic model is responsible for magnetically-induced 
P in this material. 
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Since the long magnetic modulation period [TM3 m 
Cu20Se03 means the almost collinear spin arrangement 
within a crystallographic unit cell, the experimentally ob- 
tained P-M correspondence in the collinear spin state 
(Fig. 2) should also define the relationship between local 
electric polarization and local magnetization in the SkX 
state. Figure 4 indicates the spatial distribution of local 
magnetization m{r) (Fig. 4(a)), as well as local electric 
polarization p(f) oc • m(r)) 2 e^ (Figs. 4(b)-(d)) 

and local electric charge p(f) oc V • p(f) (Figs. 4(e)- 
(g)), calculated for the SkX state given by Eq. (j3]) with 
various directions of H. The obtained p(f) and p(f) pro- 
files suggest that each skyrmion particle locally carries 
electric quadrupole moment for H || [001] (Figs. 4(b) 
and 4(e)), or electric dipole moment along the in-plane 
([001]) and out-of-plane ([111]) direction of spin vortices 
for H || [110] (Figs. 4(c) and 4(f)) and H || [111] (Figs. 
4(d) and 4(g)), respectively. These pictures are consis- 
tent with P = for H || [001], P || [001] for H || [110], 
and P || [111] for H || [111] experimentally obtained by 
the macroscopic measurement in the SkX state (Figs. 
3(c), 3(f), and 3(i)). Such a local coupling between a 
skyrmion and electric dipole strongly suggests that each 
skyrmion particle in Cu2 0Se03 can be independently 
driven by the spatial gradient of external electric field. 
The net charge is always zero within a single skyrmion 
(Figs. 4(e)- (g)), which implies the non-dissipative nature 
of skyrmion motion. 

In summary, we have investigated the magnetoelectric 
response of Cu20Se03 under various magnitudes and di- 
rections of magnetic field. The observed development of 
P in ferrimagnetic, helimagnetic, and skyrmion crystal 
spin state can be consistently explained by the d-p hy- 
bridization model, highlighting an important role of rela- 
tivists spin-orbit interaction on the magnetoelectric cou- 
pling in Cu20Se03. In the crystallized form, skyrmions 
are found to induce electric polarization along either in- 
plane or out-of-plane direction of the spin vortices de- 
pending on the applied indirection. Thus, the skyrmion 
crystal state in the insulating chiral magnet is endowed 
with the density- wave nature of both spin and polariza- 
tion. Our analysis shows that each skyrmion particle 
can locally carry electric dipole or quadrupole, which 
strongly suggests the possible manipulation of individ- 
ual skyrmion particles by external electric field (not cur- 
rent). Since electric field in insulators causes only negli- 
gible joule heat loss compared to the current-driven ap- 
proach in the metallic system, the presently established 
" magnetoelectric" skyrmion may contribute to the design 
of novel spintronic devices with high-energy efficiency. 
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